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We propose a novel way to image dynamical properties of nonstationary electron systems using 
ultrafast resonant x-ray scattering. Employing a rigorous theoretical analysis within the framework 
of quantum electrodynamics, we demonstrate that a single scattering pattern from a nonstationary 
electron system encodes the instantaneous interatomic electron current in addition to the structural 
information usually obtained by resonant x-ray scattering from stationary systems. Thus, inelas¬ 
tic contributions that are indistinguishable from elastic processes induced by a broadband probe 
pulse, instead of being a concern, serve as an advantage for time-resolved resonant x-ray scatter¬ 
ing. Thereby, we propose an approach combining elastic and inelastic resonant x-ray scattering for 
imaging dynamics of nonstationary electron systems in both real space and real time. In order to 
illustrate its power, we show how it can be applied to image the electron hole current in an ionized 
diatomic molecule. 

PACS numbers: 78.70.Ck, 42.50.Ct, 82.53.Xa, 87.15.ht 


I. INTRODUCTION 

Electron dynamics in valence shells of atoms and 
molecules determine various physical processes, such as 
chemical reactions, cooperative phenomena in solids, 
charge migration in biological systems etc. Imaging of 
electron dynamics both in real time and real space is one 
of the most important goals for modern ultrafast science 
[1-9]. X-ray free-electron lasers are a promising tool to 
achieve this task [10-14]. Wavelengths of hard x rays 
provide angstrom resolution, which corresponds to inter¬ 
atomic distances in molecular structures and solids. At 
the same time, free-electron laser sources are able to pro¬ 
duce ultrashort high intensity pulses, which give access 
to femtosecond time scales [15, 16]. 

While technological developments approach the possi¬ 
bility to make ‘electron movies’, it is necessary to answer 
fundamental questions about the interaction between a 
nonstationary electron system and an ultrashort light 
pulse. In this paper, we analyze time-resolved diffraction 
imaging by resonant hard x-ray scattering from nonsta¬ 
tionary electron systems and demonstrate that the in¬ 
stantaneous interatomic electron current is encoded in 
the Fourier transform of a scattering pattern. Resonant 
x-ray scattering (RXS) is a powerful technique that pro¬ 
vides insight into charge, orbital and spin degrees of free¬ 
dom [17-20]. RXS is an element specific method, since 
it involves transitions from atomic core shells with res¬ 
onant excitation energies that strongly depend on the 
atomic species. The resonant nature of this process al¬ 
lows one to considerably enhance the scattering cross sec¬ 
tion in comparison to the nonresonant case studied in 
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Refs. [21, 22]. This is particularly relevant for the mea¬ 
surement of valence electron dynamics in heavy elements, 
where the vast majority of electrons are stationary, since 
nonresonant x-ray scattering probes simultaneously elec¬ 
trons involved in the dynamics and electrons that are 
essentially stationary [23]. In resonant scattering, one 
can selectively enhance the scattering contribution from 
those (quasi-)particles that are actually moving. Reso¬ 
nant conditions are also essential for magnetic scattering, 
since magnetic interactions with light are very weak [24]. 
Time-resolved diffraction by resonant x-ray pulses has 
been used to reveal various ultrafast phenomena, such as 
melting of orbital and spin orders in strongly correlated 
materials [25, 26], laser-induced spin reversal [27] and 
demagnetization [28]. Applications to coherent electron 
dynamics have not yet been reported. 

X-ray Raman scattering has been proposed as a spec¬ 
troscopic probe of valence excited states in molecules 
[29, 30]. It has been suggested that information about 
valence electron dynamics can be obtained by analyzing 
Raman spectroscopy signals, such as the change in trans¬ 
mission of a probe pulse with respect to a pump pulse 
[31]. We propose a different method for measuring elec¬ 
tron dynamics, which is momentum- and time-resolved 
resonant diffraction by hard x rays. The advantage of 
this method is that one can directly image in both real 
time and real space electron dynamics such as charge 
migration, electronic wave packets in molecules, charge 
transfer during chemical reactions, as well as collective 
electron excitations. 

The differential scattering probability (DSP) of RXS 
from a stationary crystal is dominated by elastic scatter¬ 
ing processes, since their amplitudes sum up coherently 
giving rise to charge or magnetic Bragg peaks. There¬ 
fore, a most straightforward approach to calculate the 
DSP from an electron wave packet is to consider only 
the contribution of the elastic scattering processes, which 



2 


do not change the state of the wave packet. However, it 
has recently been demonstrated that inelastic (Compton- 
type) processes considerably affect scattering patterns 
obtained by ultrafast nonresonant x-ray scattering [21] 
and electron diffraction [32], leading to a loss of informa¬ 
tion about the instantaneous electron density. Thus, we 
develop a description of RXS from nonstationary elec¬ 
tron systems based on quantum electrodynamics (QED), 
since it allows taking into account both elastic and in¬ 
elastic processes correctly [33, 34]. Although the role 
of inelastic processes for time-resolved nonresonant x- 
ray scattering has been already analyzed [21, 23, 35, 36], 
these studies cannot be applied to the resonant ca.se, since 
high-energy resonant and nonresonant x-ray scattering 
are determined by different terms of the light-matter in¬ 
teraction Hamiltonian. Moreover, we find that in con¬ 
trast to nonresonant x-ray scattering [21], time-resolved 
RXS allows resolving the direction of electron flow with¬ 
out losing the connection to structural information de¬ 
spite inelastic contributions. Thereby, we show that al¬ 
though inelastic scattering processes are usually used for 
a spectroscopy analysis, and only pure elastic RXS is used 
to image structural information in stationary measure¬ 
ments, inelastic contributions serve as an advantage for 
time-resolved RXS and provide additional insights into 
electron dynamics. 

In the next Section, we derive the QED description of 
the DSP of time-resolved RXS that takes into account 
both elastic and inelastic processes and compare it to a 
“quasi-stationary” description that assumes that the con¬ 
tribution from elastic processes to a time-resolved scat¬ 
tering pattern dominates. In Section HI, we show the 
connection between the Eourier transform of a scattering 
pattern and the interatomic electron current at the time 
of measurement. We illustrate our results by describing 
a possible experiment in Section IV, showing how the in¬ 
stantaneous electron current in ionized Br 2 molecule can 
be imaged. 


the probe and the pump steps separately (see [37]). Since 
the goal of this paper is to describe how one can mea¬ 
sure electron dynamics of a given nonstationary electron 
system, the specific pump process giving rise to |4'(t)) 
is of no concern here. In Section H C, the description of 
scattering patterns from a coherent wave packet will be 
generalized for a statistical mixture of states. 


A. Quasi-stationary description 


The semiclassical approach to calculate the DSP of 
RXS from a wave packet is to substitute the wave 
packet state at the time of measurement tp, |'I'(tp)) = 
X]/C'/e“*^^‘p|4)/), for the ground state in the general 
relation for stationary RXS [38]. Then, 
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where Q is the scattering vector, c is the speed of light, 
/in(t) is the probe pulse intensity. Win is the photon en¬ 
ergy of the incoming beam, {E) is the mean energy of 
the electron wave packet and Cin is the mean polariza¬ 
tion vector of the incoming beam (we use atomic units 
for this and the following expressions). The sum over Sg 
denotes the sum over the polarization vectors of the 
scattered photons. In Eq. (2), fc{tp) is the scattering 
amplitude of atom C situated at position Rc. Jq de¬ 
notes an intermediate state with a hole in a core shell 
of atom C, Pjc is the decay width of Jc- The spatial 
distribution of the x-ray electric field is treated within 
the dipole approximation for each absorbing atom, since 
we consider a resonant process involving transitions of 
electrons from core shells, which are highly localized in 
comparison to the resonant wavelengths considered here. 


II. DIFFERENTIAL SCATTERING 
PROBABILITY FROM A NONSTATIONARY 
ELECTRON SYSTEM 

Let us consider an electron system with Hamiltonian 
Hni, which is the many-body electronic Hamiltonian in 
the absence of an x-ray field, with eigenstates [‘I’/) and 
eigenenergies Ep. We investigate scattering patterns of 
a resonant x-ray probe pulse from an electron system 
that had been excited by a pump pulse into a coherent 
superposition of the electronic eigenstates at time t = 0. 
The time evolution of the coherent superposition is given 

by 

\^{t))=Y.Cie-^^^^\^i). (I) 

I 

We assume that the pump and probe pulses do not over¬ 
lap in time. In this way, it becomes possible to describe 


B. QED description 

Let us compare the result obtained by the quasista¬ 
tionary description in Eq. (2) to the DSP derived from 
QED. If the electron system is probed by an x-ray pulse, 
the total Hamiltonian of the whole system, matter and 
light, can be written as [39] 

E = Hm T ^ ^ Wk,shi[^ ghk,5 T Rjnt? (3) 

k,s 

where ^ and ak,s are creation and annihilation opera¬ 
tors of a photon in the k, s mode of the radiation field 
with energy Wk = |k|c. iLint is the minimal coupling in¬ 
teraction Hamiltonian between the matter and the elec¬ 
tromagnetic field in Coulomb gauge 

Hint = \ J (A(r) • ^^(r), 


( 4 ) 
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where A is the vector potential of the electromagnetic 
field, p is the canonical momentum of an electron, 
and Ip are electron creation and annihilation field op¬ 
erators. We ignore the term of the interaction Hamilto¬ 
nian determined by (which dominates for high-energy 
nonresonant x-ray scattering), since it becomes negligible 
compared to Eq. (4) in the case of resonant scattering. 
The DSP is connected to the probability P(ks) of observ¬ 
ing a scattered photon with momentum kg, which differs 
from the incoming photon momenta, by 

where V is the quantization volume. P(ks) is given 
within the density matrix formalism [40] by 


P(kg) = lim Tr 

tf —^-|-oo 


( 6 ) 


where the operator 


Ok. = 5] W(wkj5]5] |<i>F;K})(<l>;^;{n'}| (7) 

S. = l F {«'} 


describes the observation of a photon in the scatter¬ 
ing mode kg, W(u;kJ represents the spectral acceptance 
range of the photon detector, {n'} is a field configuration 
that has one photon in the scattering mode kg, \^f) is 
a final electronic state vector. Pf{tf ) is the total density 
matrix of the electron system and the electromagnetic 
field at time t / after the action of the probe pulse, which 
we evaluate within the second-order time-dependent per¬ 
turbation theory using H^t as the perturbation. We take 
into account that the probe pulse duration, must be 
much shorter than the time variation of the wave packet 
during the action of the probe pulse (see Eq. (A9)). That 
means that the spectral bandwidth of the probe pulse 
must be much larger than the maximum energy split¬ 
ting among the electron states involved in the dynamics. 
Therefore, it becomes impossible to separate elastic and 
inelastic scattering events through the spectroscopy of 
the scattered photon. 

We show in Appendix A that, with those approxima¬ 
tions, the DSP of a probe pulse with intensity /in(t) = 
/q ^ which arrives at time tp after the 

pump pulse can be written as 


(W 
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where Aujj^f = Ej^—Ep and Q.f = Wk. — Win-|-AF —(A). 
Equation (8) as well as Eq. (2) describe a process where 
the system is brought into some intermediate state Jc 
by absorption of a photon, and then to some final state 


by spontaneous emission. However, the key difference 
between the expressions is that the wave-packet state 
|4'(tp)) enters Eq. (8) only once. The quasi-stationary 
description in Eq. (2) applies an analogy to an elastic 
process in a stationary system and assumes that the sys¬ 
tem is brought back into the same electronic wave-packet 
state after the scattering process. However, this analogy 
is not correct, since the initial state of a nonstationary 
system is a superposition of its eigenstates, and the elas¬ 
tic process resulting in the final state that is exactly the 
same superposition of states cannot be spectroscopically 
distinguished. This can also be shown from another per¬ 
spective. The absorption process, bringing the system 
into some intermediate state, destroys the wave packet. 
Thus, the information about the wave packet is contained 
only in the absorption term (<I)jp|ein • r|4>(tp)), but not 
in the spontaneous emission term, ($F|es How¬ 

ever, the emission step is not possible, if absorption has 
not taken place. Thereby, the emission indirectly de¬ 
pends on the wave packet and a scattered photon indeed 
provides the information about the nonstationary elec¬ 
tronic system. 

It is shown in Appendix A 2 that Equation (8) applied 
to a stationary system at Tp ^ l/Ej^ and IT(a;kJ = 1 
goes over into the conventional relation for RXS, which 
takes into account elastic and inelastic contributions [38]. 
Although Eq. (8) may resemble the expression for sta¬ 
tionary inelastic RXS, there are fundamental differences 
between them. In a time-resolved measurement, where 
the pulse duration must be much shorter than the char¬ 
acteristic time of the electron dynamics, the incoming 
photon energy is not precisely defined, and the signal 
depends on the spectrum of the broadband probe pulse 
[36]. Spectrally resolved inelastic x-ray scattering from 
a stationary system encodes a dynamic structure factor, 
which can provide access to electron dynamics at the level 
of linear response theory [41, 42]. Our method measures 
the state of an electron system at a given time indepen¬ 
dently from how this state had been created. Another 
difference is that, in contrast to the inelastic x-ray scat¬ 
tering technique [43], energy resolution is not required in 
our method in order to obtain information about electron 
dynamics, as will be shown in Section HI. 


C. General expression 


If the state prepared by the pump pulse is a statistical 
mixture rather than a coherent wave packet, the system 
at time t must be described by the density matrix 


( 9 ) 

I,K 
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In this case, the expression for the DSP becomes 

dP ^ iQ.(Rc,-RcD 

dfl 41n2c^ 

(-'q i^r 

X '^^'^IK{tp){^K\Qqr\^l), (10) 

I,K 

where Qqr is provided in Eq. (A13). This expression gen¬ 
eralizes Eq. (8) and reduces to Eq. (8) when TiK{t) = 
which is the condition that /5™(t) de¬ 
scribes a perfectly coherent wave packet. 


III. CONNECTION TO THE PROBABILITY 
CURRENT DENSITY 


In contrast to the quasi-stationary description in 
Eq. (2), which straightforwardly generalizes the elastic 
x-ray scattering theory, the DSP from a nonstationary 
system in Eq. (8) is not determined by the structure fac¬ 
tor at the time of measurement, fc{tp)P^'^^ ■ Thus, 
let us consider what information is actually encoded in 
such a scattering pattern. First, structural information 
is still present due to the factors e**^ *^*^. Second, as will 
be shown below, the probability current between scatter¬ 
ing atoms at the time of measurement is encoded in a 
scattering pattern. 

For the general situation in Eq. (9), the probability 
current density at the time of measurement, given by 

= ^Tr|p“(tp) (11) 

can be decomposed into intra-atomic and inter-atomic 
contributions as 

j(Dtp) = H ja'6“(j',^p), (12) 

Ca Ca,Cb^Ca 


where Ca and Cb refer to all atoms in the system (see 
Appendix B). The volume-integrated probability current 
between scattering atoms Cq and Cr is given by 



d^rj 


inter 

qr 


(r,tp) ■ riqr 
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= Im I (tp)(<I>j^:|(?gr|^j) 


. I.K 


where ngr is the unit vector pointing from site Cq to site 
Cr, and Gqr = J d^r|J,^(r)(V • llqr )|c’,(r), the operator 
annihilating a particle at site Rc- 
We find that the interatomic probability current is en¬ 
coded in the Fourier transform of the DSP, 


Namely, the imaginary part of the Fourier transform of 
the general expression (10) for the DSP, 

Im(J'DSp(i')) OC Jqr{tp)S[r- (Rc, - RcJ], (15) 

Cg,Cr 

has delta peaks at the interatomic distances between 
scattering atoms Cq and Cr weighted by the factor 


Jqr{tp) = Im I Y,d-IK{tp){^K\Qqr\^l) 
,I,K 


(16) 


where Gqr = dc Sdcq [see Eq. (A13) for dc,.,. and S]. 
The time-dependent factors Ijx{tp) entering in Jqr{tp) 
and jqr{tp) are weighted in both cases by matrix elements 
that are determined by the amounts of charges at sites 
Cq and Cr- Namely, the matrix element |Gqj.|<I> 7 ) 
depends on the amount of charges at sites Cq and Cr 
due to the operators ^c, and ^Cr that single out these 

charges. The operators dc, and dc,, in {^K\Gqr\^i) se¬ 
lect the charges via dipole matrix elements of transitions 
to and from inner shells localized at sites Cq and Cr that 
are determined by the charges at these sites. The op¬ 
erator S gives a time-independent contribution that can 
be factored out from the sum over / and K and, con¬ 
sequently, does not influence the temporal behavior of 
Jqr{tp). As a result, the temporal evolution of the factor 
Jqr{tp) substantially reproduces the temporal behavior of 
the interatomic current jqr{tp) (see the discussion in Ap¬ 
pendix B 1 for additional details). This relation between 
the factor Jqr and the current jqr is analogous to the rela¬ 
tion between the scattering amplitude /c and the charge 
at site C in stationary elastic RXS. In the next Section, 
we provide an example showing the connection between 
Jqr{tp) and jqr(tp) in an ionized diatomic molecule. 


IV. APPLICATION TO A DIATOMIC 
HOMONUCLEAR MOLECULE 

In order to demonstrate the basic properties of the 
time-resolved scattering patterns and compare the re¬ 
sults of Eqs. (2) and (8), we apply our formalism to RXS 
from an electronic wave packet in a diatomic homonuclear 
molecule. In recent years, significant progress in attosec- 
ond science has made it possible to prepare in a controlled 
way electronic wave packets in molecules with charac¬ 
teristic time scales of femtoseconds to sub-femtoseconds. 
It is possible to launch a wave packet by removing one 
or two valence electrons from desired orbitals [2-6] and 
even control the outcome of a simple chemical reaction 
[44] with ultrashort pump pulses. Strong-held ionization 
with subcycle optical Held transients provides attosecond 
temporal conhnement of ionization and enables a con¬ 
sequent triggering of a wave packet with a well-dehned 
phase [45, 46]. It has been demonstrated that valence- 
electron wave packets can evolve with a high degree of 
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FIG. 1. Molecular orbital diagram of the relevant states of 
BrJ. (a) The wave packet: an electron hole in the superpo¬ 
sition of the dpTTg and 4p7r„ orbitals, (b) The intermediate 
states Ja and Js: the closed valence shell and an electron 
hole in either of the two Br Is orbitals. 


coherence for much longer than 10 fs [3]. Measuring the 
induced wave packet dynamics requires sub-femtosecond 
timing synchronization between pump and probe pulses, 
which has been experimentally demonstrated in attosec- 
ond science [3, 47]. Here, we consider a coherent wave 
packet launched in a Br 2 molecule by a photoionizing 
pump pulse and imaged by an ultrafast x-ray probe pulse. 

A molecule or a solid studied by RXS should consist of 
rather heavy elements, which have a high K (or L) edge 
energy, in order to gain a sufficiently high spatial reso¬ 
lution. For instance, resolutions better than 2.5 A may 
be obtained by pulses resonant with K edges of elements 
with the atomic number Z > 22 or with L edges of ele¬ 
ments with Z > 53 [48]. Therefore, we apply our study 
to electron hole dynamics in Br 2 , which allows reaching 
a resolution of 0.9 A at the K edge of Br. 

We treat Br 2 within the molecular orbital theory, the 
LCAO approximation, and under the assumption that 
each atom contributes one atomic orbital to form a bond. 
Thus, the molecular bonding and antibonding orbitals of 
Br 2 can be expressed as [49] 

|<(.±) = (|^a)±|^b))/v/2(l±5), (17) 

with the corresponding energies 7?+ and i?_. Here, \4>a) 
and \(j)},) are basis wave functions localized at the two 
atoms of Br 2 denoted as A and B, respectively; S = 
{4>a\4'b)- The two highest occupied molecular orbitals of 
Br 2 are the and tt^ orbitals of Br 4p character (see 
Fig. la) [50]. Let us assume that the pump pulse created 
a hole initially localized at the orbital of atom A of 
Br 2 aligned in the x direction (see Fig. 2). Such electron- 
hole localization in a molecule by a photoionizing pump 
pulse is possible as has been demonstrated in Ref. [4]. 
Then, the hole starts oscillating between sites A and B 
and at time tp is a superposition of the states (17) 

|4/(tp)) = + (18) 



Probe pulse Wave packet 



FIG. 2. Schematic representation of the pump-probe scenario 
considered. The pump pulse excites a coherent superposition 
of the hole bonding and anti-bonding states of Br 2 . A reso¬ 
nant x-ray scattering pattern is taken at time tp by the probe 
pulse. 


where C± = y^(1 ± S)/2. Strong-held ionization is 
highly selective in the sense that only the most weakly 
bound orbitals can be ionized [51]. Therefore, orbitals 
that are more strongly bound than (/)_|_ and (/)_ can be 
treated as unaffected by ionization. The incident photon 
energy of the z-polarized probe pulse propagating in the 
y direction is tuned close to the K edge of Br: Win « 13.5 
keV. Thereby, intermediate states are excited with a hole 
in either of the two Br Is core orbitals, which we assume 
to be degenerate. Thus, assuming that the resonant exci¬ 
tation hlls the valence hole, only two intermediate states 
play a role, associated with the two Br Is orbitals |^j^) 
and \4>Jb) localized at A and B, respectively (see Fig. lb). 

According to Eq. (8), the DSP from wave packet (18) 
at time tp is 


dP _ loTplDoz'^ 
dH, 8 In 2c^ 


+ (ig) 

F,s^ 

+ 2Re[(Di:’j^ •e*)(Di:’j^ • e*)] sin(( 52 ,i?,r) sin(27rtp/r)} 


X Aujjp 


ry4 -I- (wk, - AujjfY 


where Rx is the interatomic distance, T = 27r/(B+ —B_), 
Boz = j d^r'^lz(j)j^ = j d?r'^lz(t)jB, = 

f d^r(j)j^^j^-^r (j)p, \(j)p) is a molecular orbital where the 
hole is situated in the final state, Fj = Tj^p = Tj^^p and 
Aujjp = Aujj^p = Aujj^p. There are two types of con¬ 
tributions to the scattering signal. The first two, time- 
independent terms in the curly braces in Eq. (19) result 
in a constant background signal. The last term depends 
on the phase of the wave packet, 2Trtp/T, and provides a 
(5a;-dependent diffraction signal. It determines the imag¬ 
inary part of the Fourier transform of Eq. (19) from Q 
space to real space, while the first two terms determine 
the real part. The time-dependent term is nonzero for 
transitions into final states for which both Tl> pj^^ and 
are nonzero, i.e. into nonlocalizable states. Thus, 
both elastic and inelastic scattering processes involving 
final-state holes in outer shells contribute to the time- 
dependent term providing information about the wave 
packet. 
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Process 

Rate (atomic units) 

2p - Is 
3p - Is 
Ap - Is 

5.2 X 10“® 

7.3 X 10“^ 

6.4 X 10'^ 


TABLE I. Spontaneous emission rates in core-excited Br cal¬ 
culated using the XATOM toolkit [52]. 


A. Experimental considerations 

The background signal is dominated by inner-shell x- 
ray emission from the 2 p —>■ Is and 3p —> Is transitions 
contributing only to the time-independent terms. It fol¬ 
lows from Table I, which shows the spontaneous emission 
rates in core-excited Br, that the ratio of non-signal to 
signal scattered photons is about 100:1. Since the pho¬ 
ton energies due to inner-shell x-ray emission are much 
lower than the incoming photon energy, the detector win¬ 
dow function should be centered at the incom¬ 

ing energy, Win, and be sufficiently narrow to suppress all 
Win in order to decrease the background. 

Another possible source of the background signal is 
due to excitation of core electrons into unoccupied states 
above the wave packet states. Therefore, the spectral 
width of the probe pulse must be small enough not to 
induce transitions into higher lying orbitals that do not 
contribute to the wave packet. But at the same time, the 
spectral width should be much larger than the energy 
splittings of the eigenstates comprising the wave packet, 
so that the probe-pulse duration is shorter than the char¬ 
acteristic time scale of dynamics in this wave packet. For 
our system, a probe pulse with a duration of 200 as still 
does not induce transitions into higher lying states of 
Br 2 and is much shorter than the period of the considered 
wave packet, T « 2 fs [50]. Scattering from molecules un¬ 
affected by the pump pulse and hence not ionized, does 
not contribute to the patterns, since there is no vacancy 
in the valence shell in this case so that resonant absorp¬ 
tion is forbidden (see Fig. 1). Only a small fraction of 
absorbed photons lead to direct photoionization of or¬ 
bitals above the K shell. Fig. 3 shows the photoioniza¬ 
tion cross section and total photoabsorption (including 
photoionization and resonant photoexcitation) cross sec¬ 
tion for Br calculated using the XATOM toolkit [52]. It 
follows from Fig. 3 that, for photon energies near the 
Br K edge, the photoionization cross section is negligible 
in comparison to the cross section for excitation of an 
electron from the Is shell to the Ap shell in Br. 

We now estimate the number of scattering patterns 
that is necessary to obtain a sufficiently high signal-to- 
noise ratio. Let us assume a peak x-ray pulse intensity 
/q = 10^® W/cm^ and pulse duration of 200 as corre¬ 
sponding to 10® photons in a focal area Aq of 1 



Incoming photon energy, eV 


FIG. 3. Photoionization cross section and total photoabsorp¬ 
tion cross section of Br calculated using the XATOM toolkit 
[52]. 


The solid angle corresponding to an independent pixel is 

dOp = = 0.15, (20) 

where A is the wavelength of the probe pulse and is the 
interatomic distance of Br 2 . We find that {dP/dVl) « 2 x 
10“® [see Fig. 4 and the discussion of it in the next Sub¬ 
section] and, thus, approximately Np = dflp x {dP/dfl) = 
3 X 10“® photons per pixel are scattered from a single 
molecule for each probe pulse. The preparation of an 
ensemble of aligned molecules provides a significant sig¬ 
nal enhancement due to the possibility to average many 
identical patterns [53]. An extremely high degree of 
molecular alignment can be achieved with a strong laser 
field (for instance, (cos^ 620 ) = 0.97 in Ref. [54] and 
(cos® 021)) = 0.89 in Ref. [55]) [56]. Alignment of Br 2 
molecules has been discussed in Refs. [57-59] and demon¬ 
strated experimentally in Ref. [60]. Recently the feasibil¬ 
ity of x-ray diffraction imaging of laser-aligned gas-phase 
molecules was demonstrated [55]. 

We assume a molecular beam density of M = 10^® 
cm“®, a molecular beam size of Im = 0.4 cm, an interac¬ 
tion area of Aq = 10“® cm® and 10% of the molecules in 
the wave packet state, so that M x x Aq x 10% « 4 
molecules per shot contribute to a scattering pattern. 
Thus, roughly 10^ shots are necessary to obtain a signal 
of 1 photon per pixel. The necessary number of shots can 
be reduced to 10® using the method in Ref. [61], which al¬ 
lows reconstructing a structure of a single molecule with 
a mean photon count on the order of 0.01 photons per 
pixel. Thus, this experiment would, in principle, be feasi¬ 
ble with the forthcoming European XFEL facility, which 
will provide 27 000 x-ray pulses per second [53]. The only 
obstacle is that the European XFEL will initially not pro¬ 
duce sub-femtosecond pulses, but this limitation can be 
overcome with the strategies described in Refs. [62-66]. 
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QED Hole density “Quasi-st.” 



- 3 - 1.5 0 1.5 3 - 3 - 1.5 0 1.5 3 

Qx, Qx, A"* 


FIG. 4. DSP from a single molecule in the Qx — Qz plane 
at Qj, = 0 according to Eq. (19) (left) and Eq. (2) (right) 
and the corresponding schematic representation of the hole 
density (middle) at different tp. Br 2 parameters: Hole lifetime 
is 260 X 10"^® s [67], Rx = 2.3 A, T = 1.7 fs [50]. Probe pulse 
parameters: Tp = 200 as, Win « 13.5 keV, 7o = 10^® W/cm^. 
A'’ is a positive integer. The dependence on polarization is 
not shown. The ranges are limited by + Q\ < 


B. DSP from the oscillating electron hole in Br 2 

Figure 4 shows the DSP from the oscillating electron 
hole in Br 2 at different times tp according to Eqs. (19) 
[QED description] and (2) [“quasistationary” descrip¬ 
tion]. Eor Eq. (19), we assumed that lE(a;kJ is centered 
at Win and its width is 5 eV, so that scattering only in¬ 
volving the Is - dpTTg and Is - 4p7r„ transitions is detected 
(see Eig. lb). This assumption influences only the ampli¬ 
tude of the DSP, but not the basic features. The signals 
are constant when the hole is localized on one of the 
atoms (at times 0 and T/2), since the scattering chan¬ 
nel for the other atom is blocked at these moments, and 
there is no interference. A diffraction signal is obtained 
when the hole is delocalized, whereby the difference of 
the positions of the maxima is equal to 2 Tr jR^. in both 
cases. The unidirectional structure of the DSP reflects 
that the molecule is aligned along the x axis. Thus, at 
times T /4 and 3T/4 in Eigure 4, the scattering patterns 
encode information about the alignment, orbital direc¬ 
tion and interatomic distances of Br 2 . 

Surprisingly, the scattering patterns derived from QED 
depend on whether the hole is moving from site A to site 


tp = r/4 tp = 3r/4 
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PIG. 5. Imaginary part of the spatial Fourier transform of the 
QED scattering patterns at times T/4 and 3r/4. The inter¬ 
ference fringes are due to the finite regions of the integrated 
functions. 


B or vice versa. They are shifted relative to each other 
by IT/Rx and are not inversion-symmetric with respect 
to Q = 0. Thus, a single scattering pattern provides 
the direction in which the electron hole is moving. In 
contrast, the DSP according to the “quasistationary” de¬ 
scription is the same at times r/4 and 3T/4, when the 
instantaneous hole density is the same. This discrepancy 
demonstrates that the “quasistationary” picture does not 
provide a correct description for a dynamical system. It 
fails because both elastic and inelastic processes are deci¬ 
sive for time-resolved measurements, for which the probe 
pulse bandwidth is much larger than energy splittings of 
electron states. 

The imaginary parts of the Fourier transforms of the 
QED scattering patterns are shown in Fig. 5 at times 
tp = T/A and tp = STjA. The amplitude and direction 
of the instantaneous interatomic electron current may be 
reconstructed by Fourier analysis of the scattering pat¬ 
terns. Let us calculate the imaginary part of Fourier 
transform of the DSP in Eq. (19): 

Im(J^DSp) «: sin(27rtp/r) {iS{x + Rx) - iS{x - Rx)) 

x^Re[(D„, .e:)(D„,.e:)]. (21) 

F 

Now, let us connect it to the electron hole current flowing 
from atom A to atom B, jabitp), and the electron hole 
current flowing from atom B to atom A, jba{tp), derived 
in Appendix B 2: 

Im(J'DSp) OC jab{tp)5{x - Rx) + jba{tp)5{x + Rx) (22) 

jabitp) = -jbaitp) = J d^rRe{4>l\Jx^a)- 

Thus, the right delta peak appears at the position cor¬ 
responding to the vector pointing from atom A to atom 
B and is proportional to the value of the electron hole 
current flowing from atom A to atom B, jabitp) (see 
Eq. (B12) for details). The left peak is proportional to 
the value of the electron hole current flowing from atom 
B to atom A, jbaitp), which is opposite to jabitp)- At 










time tp = T/4, when the electron hole is fully delocal¬ 
ized and flows from A to B [see Fig. 4], the right peak 
reaches its maximum. The peaks are switched at time 
tp = 3T/4, when the electron hole flows in the opposite 
direction. The amplitudes of the peaks are proportional 
to the amplitudes of the currents. They decrease when 
the hole localizes on one of the atoms and are zero at 
times tp = 0 and tp = T/2. 


V. CONCLUSIONS 

We have demonstrated that inelastic scattering, which 
unavoidably contributes to a scattering pattern obtained 
by an ultrashort probe pulse, can be used to one’s ad¬ 
vantage for imaging by time-resolved RXS. Applying a 
thorough theoretical analysis based on QED, we obtained 
that time-resolved scattering patterns are not determined 
by, but still connected to the instantaneous electron den¬ 
sity and contain additional information about electron 
dynamics, resolving the direction of the electron flow. 
Thereby, due to the contribution of inelastic processes, 
information resolved in a scattering pattern in a time- 
resolved measurement differs from the information ob¬ 
tained under stationary conditions. This result is quite 
counterintuitive, since one may expect that a single snap¬ 
shot from a dynamical system contains exactly the same 
amount of data as that from a stationary system. Our 
study demonstrates that such an interpretation is not 
correct, and when interpreting results obtained by time- 
resolved RXS, one should take into account that scat¬ 
tering patterns are affected by inelastic processes and do 
not coincide with those obtained in a stationary measure¬ 
ment. We have shown that time-resolved RXS provides 
additional information about the dynamical properties of 
a nonstationary electronic system along with the instan¬ 
taneous geometry, interatomic distances and directions of 
orbitals or spins within a sample. We have demonstrated 
how the direction of electron current can be retrieved 
from a scattering pattern obtained by RXS using Fourier 
analysis. 

Preparing and watching coherent electron wave pack¬ 
ets in complex matter is an emerging opportunity in at- 
tosecond science. Current probes in attosecond science 
provide access to the time domain, but do not provide 
direct access to the spatial domain. We expect that time- 
resolved RXS will prove to be a powerful tool for time- 
resolved imaging providing access to the electronic mo¬ 
tion that is the key to understanding various phenomena 
in complex molecular, biological and solid-state systems. 
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Appendix A: Details of the QED derivation of the 
DSP 

We study the resonant x-ray scattering from an elec¬ 
tron system, which is described by the density matrix 

p’^ito) = '^IiKito)\^i){^K\ (Al) 

I,K 

at time to after the action of a pump pulse. A scattering 
process consists in the annihilation of one photon and in 
the creation of one photon. Thus, the lowest-order con¬ 
tributing process demands that the operator A, which 
is linear in the creation and annihilation operators, acts 
twice. Therefore, we evaluate the total density opera¬ 
tor within the second-order time-dependent perturbation 
theory using H-mt as the perturbation: 

{n},{n},I,K 

(A2) 

Here, {n} and {n} are complete sets that specify the 
number of photons in all initially occupied field modes 
with a distribution is the second- 

order wave function at time tf, which is an entangled 
state of the electronic and photonic states: 

^/) =- f dt' (A3) 

Jto 

Jto 

where Hq = Ffm + ^k.sd^^^dk.s- to and tf are times 

prior to and after the action of the probe pulse, thus, 
to —>■ —oo and ty —)■ -l-oo. 

Probability P(wkJ according to Eq. (5) is 

2 

^(cckj = 5] IF(a;kJ 

(A4) 

= E »■("!..) E E E E (A5) 

Ss —1 {n'},F {n},{n} ki,k 2 ,si,S 2 

JC -T t \ I ^TT / 27r 27r 

xPMd"} 

ntf ntf pt2 

X / dt[ / / dt'l / d^ 

Jto Jto Jto Jto 

X |d>/) 

X (K}l4„.sXki.silM)({«}l42.S2ak,.sJ{n'}), 
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where 


1. Electric dipole approximation 


Tw,s = J ■ p)V'(r) (A6) 


are transition operators. Below, we apply that 

^'fn}.{n} (A7) 

{"}.{"} {"'} 

X ({^}|ak2.s2“ks.sj{’^'}) = Tr[/3r^aj^^_^^aki,sJ, 


where pA = E{„}.{fi} pf„},{K}l{’^})({«}| is the initial 
density operator of the radiation field [39, 40]. 

We assume that the bandwidth and the angular spread 
of the x-ray pulse are sufficiently small to satisfy the con¬ 
ditions « Win, ki ,2 ~ kin and Cki 2 ~ £in, where 

Win, kin and ein are the mean values of the photon en¬ 
ergy, wave vector and polarization of the incident beam, 
respectively. Therefore, 


27rw 


- S 

ki,k2,Sl,S2 


Tr[/5i^ai (A8) 


27r ^ , , -21n2 

= — Io[ro)e 
c 


-2 In 2 




Here, the probe pulse has a Gaussian shape and 
the amplitude of the electric field is E{ro,t) = 

. - -2 In 2 f ^ 

y'(87r/c)/o(ro)e \ J , Tq is the position of the ob¬ 

ject, tp is the time of the measurement, Tp is the pulse 
duration (FWHM of the pulse intensity) and /o(ro) = 
cE‘^{ro,t = 0)/(87r). Note that expression (A8) is the 
first-order radiation field correlation function [39, 68]. 
Since the probe pulse must be much shorter than the 
time variation of the electron density during the action 
of the probe pulse, the following condition must be sat¬ 
isfied for an appropriate time-resolved measurement 


^ EiE)(t2-ti)p^^tp), (A9) 

where p^{tp) = x^iK{ip)\^i){^K\ is the electron 

density at the time of the measurement and {E) is 
the mean energy of the nonstationary electron system. 
Therefore, the DSP at the time of measurement tp is 


(W 

dn 


f 


27rw? 


2 f-oo 

/ dwk,Wk,VP(wkJ 


(AlO) 


xEE IlK{tp){^K\f^\^F){^F\f\^l) 


F I,K 

Ef , E' „ -2 In 2 

/ dt' / dfe 

Jto Jto 


{^y 


'^k,. 






Note that |4 >f), [$/) and can be either the same or 
different eigenstates of Hm- 


We consider a resonant process involving transitions 
of core electrons to valence shells. Core shells are highly 
localized in comparison to the resonant wavelengths con¬ 
sidered here. The size of a core shell scales with the 
effective nuclear charge Z^g {Z^g for K-shells is approxi¬ 
mately equal to the nuclear charge) as n^ao/Zeg, where 
oo = 0.529 A is the Bohr radius, n is the number of the 
shell [69]. Therefore, we consider the spatial distribution 
of the electric field, within the dipole approxima¬ 
tion for each absorbing atom. Thus, the term may 
be approximated by where R is the position of the 
core shell, from which the electron is excited. This leads 
to 


= _ giEpt' (All) 

Jc 




where is an intermediate state with a hole in a core 
shell localized at Re- Since a core excited state exists for 
a very short time tj^ , its line width T =1/is taken 
into account. After the integration over time and with 
the relation p = —i[r,H^], the DSP within the dipole 
approximation becomes 


dP 

dfl 


r^Io 

4 In 2c^ 


y^,1-iK(tp) 

I,K 


(A12) 


Cg,Cr- 

g^, = dljdc„ (A13) 

dc = |$Jc)(^JclDn • f, 

Jc 


5=E / dwk„Wk^W(wkJ£^|<I> f)(<I>f|£, 




Ac^jcF(^f |< • r|^Jc) 

((wk, - ^JJJcf) + *rjc/2) 


where Q = kin — kg is the scattering vector, Up = 
Wk„ — Win + Ep — {E), Aojj^F = Ej^ — Ep and r = 
f d^rip^r) r-ip(r). Here, we applied that the transition 
energies by the absorption process can be approximated 
by Win because of the resonance condition. 

If the system is in a coherent wave packet state, when 
the matrix elements of the instantaneous electron density 
can be represented aslixitp) = and 
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the Eq. (A12) reduces to 


Appendix B: Probability current density 


(W 

dD, 


4 In 2c‘* 


/ dwk,Wk,lE(a;kJ V 

do Cl „ 


/ , ^PJJcF 
Jc 

(AM) 


{‘^fK ■ r|$j^)($j^|ein ■ r|^(tp)) 

giQ-R-c 

(^k, - Aujj^f) + iTj^/2 


2 

e 


^‘‘F^P 
4 In 2 


where |«'(tp)) = X)/ 


The field annihilation (creation) operators ■0(r)(i/)'l’(r)) 
can be expanded in terms of one-particle wave functions 

</>a(i-): 


a 

(B2) 


2. Application to a stationary system 


In this Subsection, we show that Eq. (8) converts into 
the conventional equation of the DSP for stationary sys¬ 
tems. Eor a stationary system in the state G, |4>(tp)) is 
substituted by |$g) in Eq. (A14), which leads to 


dO. 


dP JdtFit)^ ^ , , 

- y awk„Wk,Wi„IV(wkJ 

„ _i do 

E‘ 


Tr^e 


Ss = l ' 

4 In 2 


2 -^/ 7^1112 


(A15) 


E 

Jc 


AujjcF{^F\et ■ r|4>j^)(4>j^|ein ■ r|$G } 
(wk„ — 


QRc 


Here, we expressed /q via the integral f dtlin(t), where 
lin(t) = is the intensity of the incoming 


beam. 

__ (i.,kB-"in + EF-EG)^T| ^^ 

The term r^e fets /(2v7rln2) becomes 

the Dirac delta function at Tp tj^. Thus, Eq. (A12) 
for a stationary system and at Tp ^ tj^ becomes 


where Ca{c\,) annihilates (creates) a particle with a wave 
function (j)a(j)- Thus, the probability current density in 
Eq. (11) can be represented as 


j(r,tp) =Im '^IiK^{^K\c\(t)*ii{Y)Vcc<l)a{r)\^i) 

\l,K a,13 

(B3) 

If one-particle wave function (jia can be represented as a 
linear combination of functions ^^(r — Rc) centered at 
site Rc, 



=EE 7a,Ci<?^i(l'- Rc), (B4) 

C i 

then j(r,tp) can be expressed as 


j(r,tp) = 


Ca.Cb 




,I,K 


(B5) 


dF 

dfl 


f dtlin(t) dcr^* 

/T 2^ 70 ’ 


Sb = 1 


where 


(A16) 


where the operator 


^c(r) 7a,a<pi(^ - Rc)Ca 

a i 


(B6) 


da^* 

dfl 



y- (4>F|e; ■r|4>j^)($j^|ein ■r|4>G) 

^ (wk, - ^JJJcf) + Jc!^ 

(A17) 


2 


X e**^ A.ljj^f 


fA(a;kJwk,(5(wk, 


Win + Ef — Eq) 


is the conventional expression for the stationary differen¬ 
tial scattering cross section by resonant x-ray scattering 
of long pulses {Tp ^ tj^ ) taking into account both elastic 
{F = G) and inelastic {F ^ G) contributions [38]. Note 
that this relation is not correct for a stationary system, 
if Tp ~ Tj ^, and one should apply Eq. (A15) in this case. 


annihilates an electron localized at atom G. Eq. (B5) 
can be decomposed into intraatomic (Co = Cf,) and in¬ 
teratomic (Co ^ Cb) contribution. Then, the volume 
averaged electron current between scattering atoms Cq 
and Gr is 


jqr{tp) = Im(^I/K^(4>K|c]3C„|$/) ^7a,C,i7^.C,.fc 

I a,/3 i,k 

X RcJ(V • n,,.)^*(r-Rcj), (B7) 


where is a unit vector pointing from site Gq to site 
Gr. 
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1. Connection to a factor J'qr{tp) 

According to Eq. (A12), a factor Jqr{tp) = 
Im K^iK{ip){^K\Qqr\^i)^ can be represented as 


J^qriip) —Icn I ^ ' Xjp[(tp^ 

\l,K 


(B 8 ) 


X (^if l<„ • ICn • r\^i) , 


where the term 




dook^uJu^W {ujkje 4 In? 

(wk, - ^ujjfY + ry4 


(B9) 


does not depend on states / and K, and, therefore, does 
not influence the temporal behavior of Jqr{tp). Here, 
we applied that the x-ray probe pulse is resonant with 
a transition energy corresponding to the excitation of 
a core shell electron, and therefore, the contributions 
of other transitions is negligible. We also assume that 
states Jc^ and Jq^ are degenerate and differ only by the 
atom, at which the electron hole is localized (as a result, 
^^JcqF = = ^<^JF and = Tjc,. = Tj), 

since the energy splittings of core-excited states would be 
much lower than the bandwidth of an ultrashort probe 
pulse. 

Let us apply the representation of a one-particle wave 
function (pa given in Eq. (B4) to Jqr{tp)\ 

PJqripp) — dlqr ^ ^ Im ( ^ ^ X/F (4*F ICgC/^H ^ Jo.. ) 
a,a’,0,P' \l,K 

l*p,ChkloL,Caid*C^k,Cr-^Cc,l,Cq, 

) Ca,Cb,i,k 


dhbk.c. = / d^rPUr - Rcb)i<n ' r)Pcoreir - RcJ, 


(BIO) 


wave function ~ R-Ca) (analogously for 

c^pCpi). The integrals dc^b^c, for Ca 7 ^ Cq and dcbk,c,. 
for Cb 7 ^ Cr are negligible in comparison to the integrals 
dCqi,Cq and dc^kfi,., therefore, 

Jqr~Hqr ^ Ini | ^ T/F IcjjC /?/1 $ ) (Bll) 

l*0,C^klo(,Cqi 


X ( / c;^?'(/)fc(r-RcJ(e*„-r)(/)core(r-RcJ 


d^rpi{v - RcJ(ei„ • ^Wcorei^ - R-cJ 


The probe x-ray pulse must excite an electron from 
a core shell to electron states, where electron dynamics 
takes place. In this case, the terms in the sum over a, (3 
in Eq. (Bll) are nonzero for the same a and /? that enter 
interatomic electron hole current between atoms Cq and 
Cr (see Eq. (B7)), and, thus, the factors 7 ^ Crk'la^Cgi 
entering J'qr and jqr would be the same. Hence, the 
factor J'qr would be exactly proportional to jqr at any 
tp under approximation that each atom contributes one 
atomic orbital to a molecular orbital, since there is only 
one term in the sum over i and k. If it is not the 
other case, Jqr would be proportional to jqr, as long as 
for every i and k, the ratio between dCgi,Cgd’p, q and 

/ dPrpKr - RcJ(V • nqr)pi{Y - RcJ is equal. 


2. Probability current density in Br 2 

According to Eq. (13), the interatomic electron hole 
currents in the system described in Section IV are 


jab = Im 




dCg^,Cg= / C^^r«^'*(l'-R-cJ(ein-r)</>core(j'-RcJ- 


= i sin( 2 Ftp/r) J d^rRe{p*bVxpa^ 

I 


(BI2) 


jba = -sm{2Ftp/T) / d'^rRe(p*a{-Va:)pb) = -jab, 


The action of the operator cl^,Ca leads to a creation of 
an electron hole in a core shell of atom Cq with a wave 
function pcoie centered at atom Cq and annihilation of 
an electron hole in a valence shell by an electron with 


where jab{tp) is the electron hole current from atom A to 
atom B, jba(tp) is the electron hole current from atom 
B to atom A, and T = 2f/{E+ — E-). The signs of the 
gradients are opposite for jab and jba, because = 
-nsA- 


[1] F. Krausz and M. Ivanov, Rev. Mod. Phys. 81, 163 [2] O. Smirnova, Y. Mairesse, S. Patchkovskii, N. Dudovich, 

(2009). D. Villeneuve, P. Corknm, and M. Y. Ivanov, Natnre 



12 


460, 972 (2009). 

[3] E. Goulielmakis, Z.-H. Loh, A. Wirth, R. Santra, 
N. Rohringer, V. S. Yakovlev, S. Zherebtsov, T. Pfeifer, 
A. M. Azzeer, M. F. Kling, et al.. Nature 466, 739 (2010). 

[4] G. Sansone, F. Kelkensberg, J. F. Perez-Torres, 
F. Morales, M. F. Kling, W. Sin, O. Ghafur, P. Johns- 
son, M. Swoboda, E. Benedetti, et al.. Nature 465 , 763 
( 2010 ). 

[5] S. Haessler, J. Gaillat, W. Boutul, C. Giovanetti- 
Teixeira, T. Ruchon, T. Auguste, Z. Diveki, P. Breger, 
A. Maquet, B. Carre, et al.. Nature Physics 6, 200 (2010). 

[6] P. Tzallas, E. Skantzakis, L. A. A. Nikolopoulos, G. D. 
Tsakiris, and D. Charalambidis, Nature Physics 7, 781 
( 2011 ). 

[7] P. Hockett, C. Z. Bisgaard, O. J. Clarkin, and A. Stolow, 
Nature Physics 7, 612 (2011). 

[8] B. Cooper and V. Averbukh, Phys. Rev. Lett. Ill, 
083004 (2013). 

[9] J. Leeuwenburgh, B. Cooper, V. Averbukh, J. P. Maran- 
gos, and M. Ivanov, Phys. Rev. Lett. Ill, 123002 (2013). 

[10] P. Corkum and F. Krausz, Nature Physics 3, 381 (2012). 

[11] K. J. Gaffney and H. N. Chapman, Science 316, 1444 
(2007). 

[12] H. N. Chapman, A. Barty, M. J. Bogan, S. Boutet, 
M. Frank, S. P. Hau-Riege, S. Marchesini, B. W. Woods, 

S. Bajt, H. Benner, et al.. Nature Physics 2, 839 (2006). 

[13] M. J. J. Vrakking and T. Elsaesser, Nature Photonics 6, 
645 (2012). 

[14] S. R. Leone, C. W. McCurdy, J. Burgdorfer, L. S. Ceder- 
baum, Z. Chang, N. Dudovich, J. Feist, C. H. Greene, 
M. Ivanov, R. Kienberger, et ah. Nature Photonics 8, 
162 (2014). 

[15] P. Emma, R. Akre, J. Arthur, R. Bionta, C. Bostedt, 

J. Bozek, A. Brachmann, P. Bucksbaum, R. Coffee, F. J. 
Decker, et al., Nature Photonics 4, 641 (2010). 

[16] B. W. J. McNeil and N. R. Thompson, Nature Photonics 
4, 814 (2010). 

[17] J. Fink, E. Schierle, E. Weschke, and J. Geek, Reports 
on Progress in Physics 76, 056502 (2013). 

[18] T. Matsumura, H. Nakao, and Y. Murakami, Journal of 
the Physical Society of Japan 82, 021007 (2013). 

[19] V. E. Dmitrienko, K. Ishida, A. Kirfel, and E. N. Ovchin¬ 
nikova, Acta Crystallographica Section A 61, 481 (2005). 

[20] S. Lovesey, E. Balcar, K. Knight, and J. F. Rodriguez, 
Physics Reports 411, 233 (2005). 

[21] G. Dixit, O. Vendrell, and R. Santra, Proceedings of the 
National Academy of Sciences 109, 11636 (2012). 

[22] F. Zamponi, P. Rothhardt, J. Stingl, M. Woerner, and 

T. Elsaesser, Proceedings of the National Academy of 
Sciences 109, 5207 (2012). 

[23] G. Dixit and R. Santra, The Journal of Chemical Physics 
138, 134311 (2013). 

[24] J. P. Hannon, G. T. Trammell, M. Blume, and D. Gibbs, 
Phys. Rev. Lett. 61, 1245 (1988). 

[25] N. Pontius, T. Kachel, G. Schiifiler-Langeheine, W. F. 

Schlotter, M. Beye, F. Sorgenfrei, C. F. Chang, 

A. Fohlisch, W. Wurth, P. Metcalf, et ah, Applied 
Physics Letters 98, 182504 (2011). 

[26] H. Ehrke, R. 1. Tobey, S. Wall, S. A. Cavill, M. Forst, 

V. Khanna, T. Garl, N. Stojanovic, D. Prabhakaran, 

A. T. Boothroyd, et ah, Phys. Rev. Lett. 106, 217401 

( 2011 ). 

[27] C. E. Graves, A. H. Reid, T. Wang, B. Wu, S. de Jong, 

K. Vahaplar, 1. Radu, D. P. Bernstein, M. Messer- 


schmidt, L. Muller, et ah. Nature Materials 12, 293 
(2013). 

[28] B. Vodungbo, J. Gautier, G. Lambert, A. B. Sardinha, 
M. Lozano, S. Sebban, M. Ducousso, W. Boutu, K. Li, 
B. Tudu, et ah, Nature Communications 3, 999 (2012). 

[29] S. Tanaka and S. Mukamel, Phys. Rev. Lett. 89, 043001 

( 2002 ). 

[30] I. V. Schweigert and S. Mukamel, Phys. Rev. A 76, 
012504 (2007). 

[31] S. Mukamel, D. Healion, Y. Zhang, and J. D. Biggs, An¬ 
nual Review of Physical Chemistry 64, 101 (2013). 

[32] H.-C. Shao and A. F. Starace, Phys. Rev. A 88, 062711 
(2013). 

[33] S. Tanaka, V. Chernyak, and S. Mukamel, Phys. Rev. A 
63, 063405 (2001). 

[34] N. E. Henriksen and K. B. Moller, The Journal of Phys¬ 
ical Chemistry B 112, 558 (2008). 

[35] G. Dixit, J. M. Slowik, and R. Santra, Phys. Rev. Lett. 
110, 137403 (2013). 

[36] G. Dixit, J. M. Slowik, and R. Santra, Phys. Rev. A 89, 
043409 (2014). 

[37] R. Santra, V. Yakovlev, T. Pfeifer, and Z.-H. Loh, Phys. 
Rev. A 83, 033405 (2011). 

[38] Y. Ma, Phys. Rev. B 49, 5799 (1994). 

[39] R. Loudon, The quantum theory of light (Glarendon 
Press, Oxford, 1983). 

[40] L. Mandel and E. Wolf, Optical Coherence and Quantum 
Optics (Cambridge University Press, Cambridge, 1995). 

[41] P. Abbamonte, K. D. Finkelstein, M. D. Collins, and 
S. M. Gruner, Phys. Rev. Lett. 92, 237401 (2004). 

[42] J. P. Reed, B. Uchoa, Y. 1. Joe, Y. Gan, D. Gasa, E. Frad- 
kin, and P. Abbamonte, Science 330, 805 (2010). 

[43] L. J. P. Ament, M. van Veenendaal, T. P. Devereaux, 
J. P. Hill, and J. van den Brink, Rev. Mod. Phys. 83, 
705 (2011). 

[44] P. Ranitovic, G. W. Hogle, P. Riviere, A. Palacios, X.- 
M. Tong, N. Toshima, A. Gonzalez-Gastrillo, L. Martin, 
F. Martin, M. M. Murnane, et al., Proceedings of the 
National Academy of Sciences 111, 912 (2014). 

[45] A. Wirth, Th, 1. Grguras, J. Gagnon, A. Moulet, T. T. 
Luu, S. Pabst, R. Santra, Z. A. Alahmed, A. M. Azzeer, 
et ah. Science 334, 195 (2011). 

[46] M. Chini, K. Zhao, and Z. Chang, Nature Photonics 8, 
178 (2014). 

[47] A. J. Benedick, J. G. Fujimoto, and F. X. Kartner, Na¬ 
ture Photonics 6, 97 (2012). 

[48] X-ray data booklet of Lawrence Berkeley National Lab¬ 
oratory at http://xdb.lbl.gov/. 

[49] G. Ballhausen and H. Gray, Molecular orbital theory: an 
introductory lecture note and reprint volume (W.A. Ben¬ 
jamin, New York, 1965). 

[50] A. W. Potts and W. C. Price, Trans. Faraday Soc. 67, 
1242 (1971). 

[51] V. S. Popov, Physics-Uspekhi 47, 855 (2004). 

[52] S.-K. Son, L. Young, and R. Santra, Phys. Rev. A 83, 
033402 (2011). 

[53] A. Barty, J. Kiipper, and H. N. Chapman, Annual Re¬ 
view of Physical Chemistry 64, 415 (2013). 

[54] L. Holmegaard, J. H. Nielsen, 1. Nevo, H. Stapelfeldt, 
F. Filsinger, J. Kiipper, and G. Meiier, Phys. Rev. Lett. 
102, 023001 (2009). 

[55] J. Kiipper, S. Stern, L. Holmegaard, F. Filsinger, 
A. Rouzee, A. Rudenko, P. Johnsson, A. V. Martin, 
M. Adolph, A. Aquila, et ah, Phys. Rev. Lett. 112, 



13 


083002 (2014). 

[561 H. Stapelfeldt and T. Seideman, Rev. Mod. Phys. 75, 
543 (2003). 

[57] C. Buth and R. Santra, Phys. Rev. A 77, 013413 (2008). 

[58] P. J. Ho and R. Santra, Phys. Rev. A 78, 053409 (2008). 

[59] P. J. Ho, M. R. Miller, and R. Santra, The Journal of 
Chemical Physics 130 (2009). 

[60] J. B. Bertrand, H. J. Worner, P. Salieres, D. M. Vil- 
leneuve, and P. B. Corkum, Nature Physics 9 (2013). 

[61] R. Fung, V. Shneerson, D. K. Saldin, and A. Ourmazd, 
Nature Physics 5, 64 (2008). 

[62] A. A. Zholents and W. M. Fawley, Phys. Rev. Lett. 92, 
224801 (2004). 

[63] P. Emma, K. Bane, M. Cornacchia, Z. Huang, H. Schlarb, 
G. Stupakov, and D. Walz, Phys. Rev. Lett. 92, 074801 


(2004). 

[64] E. L. Saldin, E. A. Schneidmiller, and M. V. Yurkov, 
Phys. Rev. ST Accel. Beams 9, 050702 (2006). 

[65] S. Kumar, H.-S. Kang, and D.-E. Kim, Applied Sciences 
3, 251 (2013). 

[66] T. Tanaka, Phys. Rev. Lett. 110, 084801 (2013). 

[67] M. O. Krause and J. H. Oliver, Journal of Physical and 
Chemical Reference Data 8, 329 (1979). 

[68] R. J. Glauber, Phys. Rev. 130, 2529 (1963). 

[69] L. Landau and E. Lifshitz, Quantum mechanics, non- 
relativistic theory (Pergamon Press, Oxford, 1958). 



